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ABSTRACT 



Aims. We investigate the effect of the electric return currents in solar flares on the profiles of hydrogen Balmer lines. We consider the 
monoenergetic approximation for the primary beam and runaway model of the neutralizing return current. 

Methods. Propagation of the 10 keV electron beam from a coronal reconnection site is considered for the semiempirical chromosphere 
model Fl. We estimate the local number density of return current using two approximations for beam energy fluxes between 4 x 10" 
and 1 x 10 12 ergcirT 2 s . Inelastic collisions of beam and return-current electrons with hydrogen are included according to their 
energy distributions, and the hydrogen Balmer line intensities are computed using an NLTE radiative transfer approach. 
Results. In comparison to traditional NLTE models of solar flares that neglect the return-current effects, we found a significant 
increase emission in the Balmer line cores due to nonthermal excitation by return current. Contrary to the model without return 
current, the line shapes are sensitive to a beam flux. It is the result of variation in the return-current energy that is close to the 
hydrogen excitation thresholds and the density of return-current electrons. 
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1. Introduction 

The ongoing study of nonthermal excitation of the flaring chro- 
mospheric plasmas has been mainly concentrated on the ef- 
fect of particle beams coming from the coronal r e connection 
site dCanfield et all [19841: iHawlev & Fisheii [19941: iFang et ail 
119931: iKasparova & Heinzell 120021: IStepan et all 120071 an d ref- 
erences therein). Ho wever, IKarlickv & Henouxl d2002l) and 
IKarlickv et al.l d2004l) recently suggested that the role of neu- 
tralizing return currents can be as important as the role of the 
primary beam itself, both fo r intensity and linear polarization 
profiles. IKarlickv et al.l d2004l) proposed a simple model of re- 
turn current formed by the runaway electrons and compared the 
rates of atomic transitions due to collisions both with the thermal 
electrons and with the electrons of the primary beam, and due to 
collisions with the return current formed by the runaway elec- 
trons. They showed that the rates due to the return current would 
dominate the collisional processes in the atmospheric region of 
Balmer line formation. However, no calculations of theoretical 
spectral line profiles were presented. 

The aim of this paper is to take a first step towards self- 
consistent modeling of the Balmer line formation with return- 
current effects taken into account. We use a semi-empirical 
model of the flaring atmosphere as a basis for our NLTE radiative 
transfer model. Then we use a standard model for electron-beam 
deceleration due to Coulomb collisions with the ambient atmo- 
sphere and combine it with the two different physical models of 
the return-current generation. We incorporate the relevant pro- 
cesses that enter the atomic statistical equilibrium equations and 
solve them with the non-local equations of radiation transfer. At 
the end, we discuss the results and validity of our models. 



2. Electron beam and return-current propagation 

We assume an electron beam that is accelerated in a coro- 
nal reconnection site and injected into the cold chromosphere 
along the magnetic field lines. During its propag ation, the 
beam evolves under the influence of several processes dKarlickvl 
Il997l) : (a) the beam generates the return current that deceler- 
ates the beam in the return-current electric field, (b) the beam 
generates the plasma waves causing the quasi-linear relaxation 
of the beam, and (c) the beam electrons are decelerated and 
scattered due to collisions with the background plasma par- 
ticles. In the following model, we neglect the plasma wave 
processes and t he return current is taken in the form of run- 
away electrons dRowland & Vlahoslll985l: Ivan den Oord 11990: 
Norman & Smith 119781). In this form, the return-current losses 
are st rongly reduced dRowland & Vlahosl 119851 : IKarlickv et al.l 
|2004|) . Thus, only collisional losses, as described by Emslie 
( 1978), decelerate the electron beam in our case. 

Let <t) = hb^b be the particle flux of the monoenergetic beam 



of the energy 



j/2, where «b is the density of the 



beam electrons, vb their velocity, and m e th e m ass of the elec- 
tron. A ccording to Norman & Smithl (Q978) and Karlick y" et al.l 
d2004l) . a fraction of background electrons a = n^/n e forms the 
current that moves in the opposite direction in order to neutral- 
ize the electric current <?<!> associated with the primary beam. We 
use «r for the number density of the return-current electrons and 
n e for the number density of the background electrons. The neu- 
tralization condition can be expressed as 



(1) 
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Fig. 1. From upper panel: Ha, H/3, and Hy disk-center line pro- 
files for T\ = (dotted line, black), T\ = 4 x 10 11 (thick solid 
line, blue), 6 x 10 11 (thick dash-dotted line, green), 8 x 10 11 
(thin solid line, red), and 1.0 x 10 12 (thin dash-dotted line, vi- 
olet) ergcm^s -1 . No return-current effects are taken into ac- 
count. Note that nonthermal profiles almost overlap in this range 
of beam fluxes. 



We assume that all the return-current electrons move with the 
same superthermal velocity vr. Note that the return-current flux 
only depends on the total flux of the beam. For a realistic power- 
law distribution of the beam, the main part of the beam flux is 
given by electrons with energy close to the low-energy cutoff of 
this distribution. For these reasons and for simplicity, we con- 
sider the monoenergetic beam in our model. 

We used two mode ls for estimating «r. First, following 
iNorman & Smith! {1978), the number of runaway electrons can 
be estimated as 



«r 1 
— = - exp 

He. 2 




(2) 



where £/£d = n B v B /n e v Tc (see iKarlickv et al.ll2004t Ixii et alj 
120051 for definition of the electric field £ generated by the elec- 
tron beam and the Dreicer electric field So), and Vy e stands for 
the thermal velocity of the background electrons. In the second 
model we assume that the return current is formed by a fixed rel- 
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Fig. 2. Same as Fig. Q] plus including collisions with return- 
current electrons. The value of a is given by Eq. (0. 



ative number of background electrons everywhere in the upper 
chromosphere: 



«R 

a — — = const. 

«e 



(3) 



In this model, a is obtained by averaging the values from the pre- 
vious model over the Balmer line formation layers. The resulting 
value for each beam flux can be found in TableQ] 

The normalized energy distribution of electrons can be lo- 
cally expressed in the form 

f(E) dE = [c M f M (E) + c B 6(E - E B ) + c R 6(E - E R )] dE , (4) 

where £r = m e v R /2 is the return-current energy and c, stands for 
the normalization coefficients (i.e., cr = «r/(«m + »b + »r), etc.). 
The index M stands for the background electrons (without the 
runaway ones), and they obey the Maxwell-Boltzmann energy 
distribution. 



3. Hydrogen-electron collisions 

To take the effect of the beam/return current into account, we 
have to calculate all the electron-hydrogen excitation rates, the 
rates of ionization by electron impacts, and the rates of the 
inverse processes. We use the data for total collisional cross- 
sections for the bound-bound and bound-free transitions by 
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Fig. 3. Same as Fig. Q]plus including collisions with return- 
current electrons. The value of a is constant along whole beam 
trajectory. 

iJanev & Smith! (Il993l) . retrieved through the GENIE database 
(http : //www-amdis . iaea . org/GENIE/). We do not consider 
any atomic polarization or angular dependence of the collisional 
processes in this work. The excitation or deexcitation rate of the 
n — > m transition between the two shells can be calculated using 
the common formula 

f°° [2E 

C nm = n e dE J — f(E)o- nm (E), (5) 
Jo V m e 

where cr nm (E) is the total cross-section of n — > m at impact en- 
ergy E, and f(E) is electron energy distribution (|4j. The cross- 
sections for deexcitation m — > n (m > n) were calculated using 
the n — > m cross se ction and the principle of detailed balance 
(e.g.. lJefferieslll968l) . 

^ffln(^) — T- 0~nm (E + E fWl ) , (6) 

8m E 

where E nm is the excitation threshold of the transition and g n / m 
the statistical weight of the given level. 

The rates of the inverse process of ionization by an elec- 
tron impact, the three-body recombination e + e + p— > Hi + e, 
must be treated separately due to the nonthermal natu re of 
the problem. Using the arguments of detailed balance dFowleri 
fl955tlJefferiel T968) in the quasi-classical approximation of the 



electron-hydrogen collisions, one can derive a formula for the 
recombination rate. Let us consider the ionization of level n 
with the ionization energy E„ by an electron with energy 
Once cross-section S nc (Ej,E a ) = dcr„ c /dE a of the encounter 
after which we find the two electrons with energies E a and 
Et, = Ej - E„ - E a is known, one can write the total recom- 
bination rate as 

C cn = Lf dE a [ dE b fiEa r ^ b) E i S„ c (E i ,E a ), (7) 
Jo Jo yE a Eb 

where L = nln 2 h 3 /167rm 2 x 6.97 x 10~ 27 n 2 n 2 in cgs units, 
and h stands for the Planck constant. The expression (0 can 
be reduced to the well-known expression Cj„ hermal = 2.06 x 
10" 16 n 2 r- 3/2 n e exp(£ , „/^ B r)Cj e hermal in the particular case of 
Maxwell-Boltzmann distribution, where T e is the temperature of 
electrons and k% the Boltzmann const ant. For the cro ss section 
S„ c we used the approximate data of Omidvar (1965) at lower 
energies and the classical formula of iThomsonl (|1912|) at high 
energies. In our case, the electron energy distribution function 
given by © leads to 9 terms that contribute to the recombination 
rates. We verified numerically that, in all our models, the rates 
of three-body recombination involving the nonthermal electrons 
of the return current are more than one order of magnitude be- 
low the rates of the same processes involving the thermal elec- 
trons. Our tests show that neglecting the nonthermal three-body 
recombination affects the resulting profiles in a very negligible 
way. The three-body recombinations involving the beam elec- 
trons are completely negligible since their rates are several or- 
ders of magnitude below the thermal ones. Thus, we took only 
the thermal term containing m 2 fy[{Ea)fy[{Et,) into account. 

4. Results 

We calculated the NLTE radiative transfer for a 5 -level plus con- 
tinuum hydrog^rnismg_tiiesMniempirical ID plane-parallel flare 
model Fl (Machad o et alJlT980h in which the temperature struc- 
ture was kept fixed; in this way, we foun d the differential ef- 
fects on the Ha, HyS, and Hy lines (cf. iKasparova & Heinzell 
2002). We used the preconditioned equations of statistical equi- 
librium dRvbicki & Humrrierll 19911) and solved the coupled sys- 
tem of NLTE equations by the ac celerated lambd a iteration 
(ALI) method. For further details, see lHeinzel (fl9 95). We found 
the equilibrium state for several beam fluxes with or without the 
return current. The initial beam fluxes chosen in our calculations 
were:^i =4xlO n , 6x10", 8x 10 11 , and 1 x lO^ergcirrV 1 . 
If the fluxes were lower, the number of runaway electrons would 
decrease fast and their energy would exceed the beam energy in 
most depths. This gives a limit for using of this simple model. 
Higher fluxes, on the other hand, would be unrealistic. The initial 
energy of beam electrons was set to £o = 10 keV. 

In Fig.Q] there are first three Balmer line profiles that result 
from the nonthermal bombardment by the primary beam. The 
effects of return current were completely ign ored. In this sense , 
these calculations are similar to the ones of lFang et all (119931) 
andlKasparova & Heinzell d2002l) . In spite of their probably lim- 
ited physical relevance, these profiles are useful for demonstrat- 
ing the effects of return currents in the more appropriate models 
that follow. Figure [2] shows the situation where a is calculated 
using Eq. {2); i.e., the relative number of runaway return-current 
electrons is calculated at each depth in the atmosphere. Finally, 
in Fig. [3] there are profiles for the model with a constant along 
the atmosphere. In the layers of Balmer line formation, a re- 
mains approximately constant and its mean values are shown in 
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Table 1. The properties of the return currents. T\ stands for the 
initial flux of the 10 keV beam, a is the mean relative number of 
runaway return-current electrons in the Balmer lines formation 
region, and Er stands for the typical energy of the return-current 
electrons in these layers. Both a and £r are average quantities 
which can roughly characterize the return-current properties in 
the region of interest. 



Ti [10" erg cm- 2 s" 1 ] 


a 


£ R [eV] 


4 


0.01 


100 


6 


0.05 


14 


8 


0.11 


6 


10 


0.16 


3.8 



TableQ] Comparing Fig.[TJwith Figs.|2]and[3] one can see that the 
effect of return current is very significant: All three lines show a 
prominent increase emission in the line center. 

In the region of Balmer line formation, the energy of the re- 
turn current can be close to the excitation threshold of n > 3 
levels of hydrogen, and it remains approximately constant along 
an extended trajectory. The beam is finally stopped on a very 
short path. The overall path of the beam is, however, sensitive 
to the initial energy of the beam. In order to model the beam 
propagation and line formati on accurately, one has to interpo- 
late the orig inal Fl model of MachadoetalJ dl980h by a number 
of grid points in the layers of the Balmer line formation. Since 
the return-current energy and density are sensitive to the beam 
flux (see Table [TJ, the resulting variation of the nonthermal col- 
lisional rates leads to a significant variation in line profiles. For 
both models under consideration, a maximum emission is found 
for the beam flux of 6 x 10 11 ergcirT 2 s , although the result- 
ing profiles from these models differ slightly from each other. 
In contrast to the case of 4 x 10 11 ergcirT 2 s -1 , for which we 
found the least emission among the studied flux intervals, the 
return-current density is higher by a factor of 5. It leads to a sig- 
nificant increase in nonthermal excitation rates. The disagree- 
ment of the profiles at fluxes below 6 x 10 11 ergcnT 2 s~' is the 
result of the significant dependency of a on the beam flux and 
atmospheric depth. The values of a at low fluxes (shown in the 
Table [TJ are only a rough approximation for the Balmer line for- 
mation layers, and the model with a = const, seems to be less 
accurate than the one given by Eq. ©. On the other hand, good 
correspondence between the models is found for high fluxes. In 
this case, the variation in a is less sensitive to the beam flux 
and does not strongly vary with atmospheric depth. Then, the 
a = const, model seems to give the appropriate results. The rea- 
son the higher beam fluxes lead to a lower emission in the lines 
is that the energy of the return current is not sufficient to excite 
hydrogen atoms as can be seen in Table Q] 

5. Conclusions and outlook 

In this paper we used a simple model of the 10 keV electron 
beam propagating in the chromosphere. We used two different 
models of the return-current formation and calculated the differ- 
ential effect on the profiles of the hydrogen Ha, H/?, and Hy lines 
of the semiempirical Fl model. 

The return-current flux only depends on the total flux of the 
beam. For a realistic power-law distribution of the beam, the 
main part of the beam flux is given by electrons with energy 
close to the low-energy cutoff of this distribution. Therefore, 
for simplicity we used the monoenergetic beam in our model. 
Moreover, the excitation and ionization cross-sections of low- 



energy return-current electrons are larger than those for beam 
electrons, which makes the return-current effects on line core 
formation stronger. Taking high-energy beam electrons into ac- 
count (i.e. using power-law distribution) would lead to increased 
emission in the line wings due to penetration of those electrons 
into the deeper atmospheric layers. However, the total flux of the 
beam in these layers is significantly lower than the initial beam 
flux, and the return current and corresponding effects are also 
strongly reduced. 

Even our simple model shows that the effect of return cur- 
rent is very important for future study of the hydrogen lines for- 
mation since the energy of the return current can be expected 
to be on the order of the excitation threshold energies of up- 
per hydrogen levels, for which the excitation cross-sections are 
high. Moreover, the fluxes O are high enough to excit e a suf - 
ficient number of atoms. As shown by Karlic kTet alj (|2004), 
the collisional rates from the nonthermal collisions can domi- 
nate the collisional rates in the Balmer line formation regions. 
The two models used in this work lead to similar results for 
higher energy fluxes, but the result differs for lower fluxes. The 
excitation threshold effects seem to play an important role for 
higher fluxes, but they are very likely only a consequence of the 
monoenergetic model we used. The difference between the two 
models shows the used approximations to be incompatible for 
lower fluxes, where the the approximation of constant a is not 
applicable. A detailed description of the energy distribution of 
the return-current electrons would lead to more realistic line in- 
tensities. This complex issue will be subject of a forthcoming pa- 
per, which will also study the impact polarization of the Balmer 
lines. 
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